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I.  INTRODUCTION 


1.1  Program  Objectives 

The  ultimate  goal  of  the  work  reported  herein  is  the  development  of 
a real-time  method  for  the  short-range  (0-20  minute)  forecast  of  storm 
development  and  motion.  The  initial  stop  in  this  program  was  the  develop- 
ment of  objective  techniques  for  the  efficient  representation  of  informa- 
tion obtained  from  a single  Doppler  weather  radar  (Crane,  1077).  The 
radar  data  are  processed  to  obtain  fixed  and  peak  referenced  reflectivity 
contours  and  peak  referenced  tangential  shear  contours.  The  essential 
information  contained  in  the  contour  data  are  represented  In'  a set  of 
attributes.  In  this  report,  we  consider  the  second  step  in  the  program, 
the  assembly  of  data  obtained  from  a set  of  radar  scans  into  a set  of 
attributes  that  characterize  the  three-dimensional  structure  of  the 
small  (tree ipi tat  ion  cells  and  surrounding  echo  regions,  their  motion,  and 
their  development  in  time.  Initial  consideration  is  given  to  the  problem 
of  forecasting  precipitation  cell  development  and  motion. 

1.2  Summary  of  Results 

The  work  performed  under  this  contract  was  concentrated  in  three 
major  areas:  (a)  improvement  to  the  cell  detection  program  developed 
under  the  previous  contract  (Crane,  1977),  (b)  development  of  the  volume 
cell  detection  and  tracking  program,  and  (c)  preparatory  data  analysis  to 
support  the  development  of  cell  forecast  algorithms. 

lhe  cell  detection  program  was  extensively  revised  to  reduce  the 
amount  of  core  storage  required  for  some  of  the  integer  arrays  used  in 
peak  referenced  contour  analysis;  to  process  reflectivity,  tangential 
shear,  and  Doppler  spread  data  simultaneously;  to  include  spread  information 
in  the  reflectivity  and  tangential  shear  peak  attributes;  to  use  the 
jelativel\  spikv  nature  of  ground  clutter  as  at:  aid  in  rejecting  ground 
clutter;  to  resolve  autumn t i on  1 1 > Doppler  velocity  ambiguities  using 
radial  velocit>  data  from  resolution  elements  adjacent  in  range  and  in 
azimuth;  and  to  van  the  quantization  steps  used  in  the  peak  referenced 
contouring  operation.  In  the  process  of  revising  the  program  additional 
refinements  were  made  in  the  cell  detection  logic  to  correct  several 


defects  detected  using  actual  radar  data  and  to  improve  the  technique 
used  to  store  temporarily  the  peak  referenced  contour  data  to  increase 
the  speed  of  the  processing  algorithm  and  reduce  overall  core  storage. 

At  present,  the  entire  program  required  for  the  simultaneous  processing 
of  reflectivity,  tangential  shear,  and  Doppler  spread  data  fits  into  the 
core  storage  available  on  the  C1K'  ooOO;  a result  not  anticipated  at  the 
start  of  the  program. 

The  program  is  currently  relatively  slow  in  operation  due  to  (11  the 
extensive  number  of  input  and  output  tape  or  disk  operations  used  to 
read  the  input  radar  data,  store  data  for  plotting,  store  data  for  sub- 
sequent volume  cell  detection  and  tracking  and  to  store  data  for  B scan 
printouts  and  (J)  the  large  number  of  bit  manipulation  operations 
required  to  pack  and  unpack  data  for  storage.  A considerable  time 
savings  can  still  be  accomplished  by  reducing  the  amount  of  output  data 
and  by  rewriting  the  program  in  assembler  language  to  reduce  the  time 
spent  in  subroutine  calls  for  bit  manipulation  and  indirect  storage 
addressing . 

A volume  cell  detection  and  tracking  algorithm  was  developed  to  com- 
bine  the  data  from  successive  radar  scans . The  fixed  and  peak  referenced 
contouring  operations  are  performed  independently  for  each  azimuth  scan 
of  the  radar.  The  reflectivity  cells  (peak  referenced!  are  associated 
from  one  elevation  angle  to  the  next  to  form  a volume  cell  and  from  one 

' 

volume  cell  to  the  next  to  form  a track.  The  reflectivity  cells  are 
used  to  identify  the  basic  organization  of  the  radar  data.  Fixed  contours 
that  enclose  associated  reflectivity  cells  are  used  to  determine  the  o- 

. | 

dimensional  echo  envelope  for  the  reflectivity  cell  (or  cells').  Tangential 
shear  and  Doppler  spread  cells  are  likewise  associated  from  scan-to-sean; 
associated  with  reflectivity  cells  (if  possible-);  and  identified  with 
fixed  contour  envelopes.  The  lowest  threshold  fixed  contour  that  encloses 
each  of  the  reflectivity,  shear,  or  spread  cells  is  used  to  tag  that  cell 
to  aid  in  the  association  processing. 

A new  set  of  attributes  are  developed  for  each  volume  cell.  A list  • 

of  these  attributes  are  given  in  Table  I for  fixed  contours  and  for  re-  1 

flectivity,  shear,  and  spread  cells.  The  attributes  to  be  calculated  for 
each  volume  cell  are  readily  modified.  The  basic  volume  cell  detection 
and  tracking  algorithm  is  used  to  perform  the  scan-to-scan  association. 
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When  data  from  two  scans  arc  associated,  they  can  he  used  readily  to 
calculate  any  desired  attribute.  The  lists  given  in  Table  I are  pre- 
liminary in  nature  covering  the  parameters  we  currently  expect  to  be  most 
important  in  subsequent  analyses,  Experience  vs  i t h the  use  of'  those  pro- 
grams for  a larger  data  set  will  be  required  before  a final  set  of  attri- 
butes can  be  specified. 

Experience  with  a larger  data  set  is  required  as  a prerequisite  to 
the  development  of  cell  forecast  algorithms.  Crane  (lO’o)  processed 
reflectivity  data  from  a number  of  radar  scans  obtained  at  Wallops 
Island  during  the  summer  of  19".>  bv  the  Johns  Hopkins  Applied  Physics 
Laboratory  (API.)  in  the  process  of  initially  establishing  the  utility  of 
the  peak  referenced  contouring  technique.  These  data  were  used  by  Crane 
ll9"o)  to  investigate  cell  lifetime  and  possible  cell  t rack ing/evt rape  1 a 
t ion  techniques.  In  this  study,  we  used  the  same  data  set  to  determine 
the  relative  spacings  between  cells.  The  object  of  the  spacing  study 
was  to  determine  if  preferred  spacings  are  evident  in  nature.  If  so, 
the  preferred  structure  can  be  used  to  forecast  the  most  probable  loca- 
tion for  new  cell  development.  The  ability  to  forecast  locations  for 
new  cell  development  would  add  significantly  to  the  0-20  minute  forecast 
because  the  median  cell  lifetime  is  less  than  JO  minutes  and  forecast- 
ing by  extrapolation  along  a cell  track  is  not  adequate.  \n  examination 
of  the  Wallops  Island  data  revealed  that  the  nearest  neighbor  distances 
between  cells  were  between  7 and  9 km  and  that  these  distance  values  did 
not  depend  on  the  type  of  rain  observed.  This  suggests  that  new  cell 
site  forecasts  are  feasible. 

1 .5  Software  Development 

The  goal  of  this  contract  with  the  Air  Force  Geophysics  Laboratorx 
(AFGL)  is  to  provide  efficient  computer  software  to  obtain  parameters  to 
represent  the  essential  information  obtained  in  a sequence  of  scans  of 
a single  Doppler  weather  radar.  The  radar  used  to  provide  the  data  is 
the  C-hand  Doppler  radar  operated  by  the  Weather  Radar  Branch  of  AFGl,  at 
Sudbury,  Massachusetts.  The  computer  programs  were  prepared  for  the 
GDC -0000  at  AFGL. 

There  are  now  three  separate  programs  in  the  sequence  of  programs 
to  be  used  for  weather  radar  data  processing.  The  first  is  the  cell 


8 


detect  ion  program  which  was  developed  under  the  previous  contract  (it. mo. 
IP"'')  and  extensively  modified  during  the  course  of  this  contract.  The 
second  is  a plotting  program  used  tv'  display  contours  and  cell  centroids. 
Vhis  program  was  developed  under  the  previous  contract  and  modit'ied  under 
this  contract  tv'  displax  separated  the  locations  of  the  centroids  of  the 
fixed  echo  regions  and  the  centroids  for  each  of  the  peak  detected  cells 
ref l ect i v i t v , tangential  shear,  and  Peppier  spread  cells.  Ihe  third  pro 
gram  detects  and  tracks  volume  cells.  It  was  developed  under  this  con 
tract.  This  program  uses  data  generated  In  the  cell  detection  program 
to  develop  three-dimensional  cells  and  to  describe  their  location,  height 
and  intensity.  The  volume  cell  attributes  are  listed  in  table  I. 

The  ultimate  goal  of  this  work  n i l 1 be  the  development  of  real  time 
radar  processing  techniques  tv'  be  used  on  a dedicated  computer  svstom 
that  is  an  integral  part  of  the  radar  svstom.  The  programs  generated  tv' 
date  on  the  current  contract  are  exploratorx  in  nature.  They  were  design 
tv'  fit  within  the  core  storage  limitations  of  the  I’lV-onOP  computer  but 
still  have  a high  degree  of  fle.vihi  l i ti  in  modifying  storage  arrav  si.es 
and  in  providing  auxilliarv  output  for  testing  the  program.  The  program 
design  had  real-time  processing  in  mind  and.  in  the  end,  should  be  reason 
able  efficient  when  tailored  tv'  an  on-site  computer  svstom.  listing  and 
operating  instructions  for  the  programs  developed  on  this  contract  arc 
included  in  the  appendices. 

1.4  Organisation  of  the  Keporf 

\ summarx  of  the  modi  fic.it  ions  tv'  the  cell  detection  program  and 
plotting  program  is  given  in  Section  .'.  \ description  of  the  volume  cell 

detection  program  is  given  in  Section  Initial  cons  iderat  ion  of  the 
cell  forecasting  problem  is  given  in  Sect  ion  I.  Sample  outputs  are  pro 
vided  in  Section  Section  o summari  cs  progress  tv'  date  and  provides 
rccommend.it  ions . 


IMPROVEMENTS  TO  THE  CELL  DETECTION  PROCRAM 


2.1  Cell  Detection  Logic  Changes 

The  cell  detection  program  was  modified  to  include  the  detection  of 
Doppler  spread  peaks  and  to  pack  integer  addresses  used  in  the  peak 
detection  routine  six  to  a CDC  word.  The  latter  change  was  required  to 
make  room  for  the  additional  data  used  in  Doppler  spread  processing.  In 
the  process  of  revising  the  address  storage  procedure,  a number  of  minor 
logic  errors  were  detected.  These  errors  were  associated  with  the  opera- 
tion of  the  radial -to-radial  association  of  contours  at  different  thres- 
hold levels  (see  Section  5.2  of  Crane,  1977).  The  errors  have  been  cor- 
rected and,  in  the  process,  changes  were  also  made  to  streamline  the 
storage  of  the  temporary  contour  data. 

At  present,  the  temporary  attribute  storage  array  (TATR  in  subrou- 
tine PF.AKD,  see  Appendix  D)  is  used  to  identify  active  cells  at  each  pos- 
sible threshold  level  (relative  to  the  peak  level  for  the  cell)  and  to 
store  pointers  to  previously  associated  active  cells  as  well  as  provide 
partial  attribute  summations  for  the  active  cells.  The  TATR  array  is 
doubly  subscripted  maintaining  temporary  storage  for  the  number  identifi- 
cation of  the  peak  for  the  active  cells.  For  each  active  cell  the  stored 
data  includes  the  current  estimate  of  the  peak  value;  a set  of  partially 
summed  attribute  values  for  each  of  the  nested  contours  for  each  possible 
threshold  level  below  the  peak,  area,  average  value,  centroid  location 
plus,  for  reflectivity  cells,  Doppler  spread,  radial  shear,  tangential 
shear,  and  radial  velocity;  the  azimuth  count  for  the  last  azimuth  on 
which  the  attributes  for  each  contour  were  updated;  and  a pointer  to  the 
enclosing  fixed  contour.  This  last  item  is  included  for  later  use  in  the 
volume  cell  detection  program.  If  a cell  is  not  active  at  a particular 
threshold  level  but  was  previously  active  on  the  current  radial,  the  area 
attribute  contains  a pointer  to  the  cell,  now  active  at  that  level.  After 
processing  the  data  for  a radial  the  cell-to-cell  pointers  are  zeroed. 

If  a cell  becomes  inactive,  all  pointers  to  that  cell  are  also  zeroed 
before  processing  data  for  the  next  radial. 

'he  program  as  it  now  sands  is  configured  to  allow  rapid  changes 
in  the  dimensions  of  all  the  storage  arrays.  This  convenience  costs  in 
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producing  relativelv  longer  processing  times  than  t'or  a program  with 
lived  storage  allocations.  \ final  program  modification  should  be  made 
when  the  exploratory  runs  of  the  computer  program  are  finished  to  fix  the 
array  sizes  and  optimize  the  program  for  rapid  operation. 

2.2  Refinements  in  Ve’ocitv  Data  Processing 

The  program  was  modified  to  provide  automatic  resolution  of  Doppler 
velocity  ambiguities.  the  first  detected  velocity  data  for  each  scan  is 
assumed  to  be  in  the  velocity  interval  spanning  zero  velocity.  Successive 
observations  at  different  but  adjacent  range  or  azimuth  locations  are 
assumed  to  differ  from  the  first  observation  by  less  than  the  velocity 
ambiguity.  If  the  difference  is  larger,  it  is  adjusted  so  that  it  is 
smaller  by  adding  or  subtracting  a value  equal  to  the  velocity  ambiguit> 
(X/2*PRF  where  \ is  the  wavelength  and  PRF  is  the  pulse  repetition  fre- 
quency!. The  data  for  each  radial  are  then  reexamined  to  ensure  that 
each  velocity  observation  differs  from  the  average  of  the  velocity  values 
for  the  two  surrounding  range  intervals  by  less  than  the  velocity  ambi- 
guity and  at.!jusfe<!  again  if  required. 

The  above  processing  is  used  for  contiguous  range  and  azimuth  ele- 
ments. If  large  spaeings  exist  between  regions  of  radar  data  above  the 
processing  threshold,  the  first  data  point  in  one  region  is  compared  with 
the  average  velocity  for  the  previous  region.  The  data  therefore  are 
objectively  adjusted  to  remove  ambiguities.  Provision  is  provided  to 
reject  data  from  this  analysis  if  the  Doppler  spread  is  not  within  pre- 
set bounds  but  this  option  has  not  been  exercised. 

The  processed  data  appear  to  be  reasonably  smoothly  varying  as 
indicated  in  Figure  1.  The  tangential  shear  data  obtained  by  numerically 
differentiating  the  radial  velocity  data  are  however,  more  variable  as 
illustrated  in  the  figure.  The  apparent  noise  evident  in  the  tangential 
shear  field  led  to  the  inclusion  of  Doppler  spread  data  in  the  analysis 
process.  The  program  was  also  modified  to  introduce  a new  constant  to 
specif)  the  quantization  interval  for  peak  detection  contouring.  By  set- 
ting this  constant,  additional  control  is  exercised  over  the  peak  detec- 
tion process.  Referring  to  Figure  l,  no  peaks  would  be  detected  if  the 
quantization  step  is  1 m/s  but  a number  of  peaks  are  detected  for  a 0.1 
m/s  quantization  step  (see  Section  51. 
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Figure  1 Saaple  Data  Set  Including  Two  Fixed  Contour  Echo  Areas  and  Three  Reflectivity  Peaks 


Doppler  Spread  Processing 


Velocity  variance  data  may  now  bo  routinely  processed  using  the  com- 
puter program,  lac  it  it i es  to  analyze  this  data  field  were  added  because 
of  the  apparently  high  noise  levels  in  the  tangential  shear  data  and  the 
expectation  that  at  high  signa 1 -to-noi se  ratios,  the  spread  data  may  be 
more  reliable.  As  discussed  before  (Trane,  h>77),  the  spread  data  arc 
assumed  to  contain  the  same  information  as  the  tangential  shear  data  due 
to  the  dominant  effect  of  velocity  fluctuations  on  the  scale  of  the 
antenna  beamwidth  at  the  observation  range.  Both  data  fields  are  now 
processed  to  allow  an  i nt ereompar i son  between  the  information  available 
from  each  field  and  to  investigate  their  relative  sensitivity  to  noise. 

Ihe  velocity  variance  output  from  the  pulse  pair  processor  is  known 
to  be  adversely  affected  by  low  signal-to-noise  ratios.  A noise  level 
dependent  threshold  parameter  is  provided  for  thresholding  the  velocity 
variance  data  but  has  not  been  exercised  as  yet.  To  reduce  the  sensi- 
tivity of  the  variance  data  vine  to  noise,  the  square  root  of  the  variance 
is  used  for  analysis  of  the  Poppler  spread.  A sample  of  this  data  is 
depicted  in  figure  1.  It  is  noted  that  the  velocity  data,  radial  velocit>, 
Poppler  spread,  and  tangential  shear,  are  processed  only  when  the  reflec- 
tivity values  are  above  the  lowest  (processing!  threshold  level  which, 
for  this  figure,  is  JO  dllZ.  High  Poppler  spreads  are  evident  in  two 
regions  of  the  figure,  (1)  at  low  reflectivity  values  and  (J1  at  high 
reflectivity  values.  The  high  spreads  at  low  values  are  presumably  caused 
by  signal-to-noise  problems  and  should  be  suppressed  using  the  signal-to- 
noise  threshold  level.  The  high  values  corresponding  to  the  higher 
reflectivity  values  are  presumably  real,  although  the  correlation  with 
the  tangential  shear  values  is  not  good  for  the  data  in  this  figure.  A 
larger  data  sample  must  be  analyzed  to  optimize  the  processing  of  tan- 
gential shear  and/or  Poppler  spread  data. 

Mod i f i ca t ions  for  Poppler  spread  processing  were  made  in  several 
areas  of  the  program.  I'he  largest  change  included  the  addition  of  tem- 
poral’' attribute  storage  arrays  for  the  detection  of  Poppler  spread  peaks. 
In  addition,  the  reflectivity  peak  and  tangential  shear  peak  detection 
algorithms  were  changed  to  increase  the  number  of  attributes  so  the  aver 
age  Poppler  spread  within  a detected  cell  could  be  determined.  This  attri 
bute  i s now  listed  with  the  others  as  described  in  Section  S. 


3.  VOl.lIMP.  Cl-LL  m-TI.CTION  AND  TRACKING 


3.1  Definition  of  a Volume  Cell 

The  cell  detection  program  provides,  as  one  of  its  outputs,  a list 
of  attributes  for  fixed  contours,  reflectivity  peaks,  tangential  shear 
peaks,  and  Doppler  spread  peaks.  At  a minimum,  each  attribute  iist 
defines  the  average  value,  area,  and  centroid  location  of  each  cell  for 
each  azimuth  scan  of  an  elevation  scan  sequence.  The  expected  scan 
sequence  for  the  radar  includes  azimuth  scans  (or  azimuth  sectors  in  a 
raster  scan)  at  successively  higher  elevation  angles  until  a complete 
volume  scan  is  completed.  The  data  for  each  elevation  or  tilt  of  the  scan 
sequence  must  be  combined  with  the  data  from  other  tilts  to  provide  a 
three-dimensional  picture  of  each  cell  or  echo  volume. 

The  process  of  combining  data  from  scan  to  scan  in  a volume  scan 
sequence  is  illustrated  schematically  in  figure  This  figure  repre- 
sents a height  section  through  three  cells.  The  contours  represent 
quantized  data  (e.g.  reflectivity  data  in  1 dB  stops)  that  have  been  pro- 
cessed on  each  azimuth  scan  to  produce  a set  of  reflectivity  peak  cells. 
The  minimum  attribute  set  is  identified  by  the  length  of  the  heavy  bar 
(detected  length  or  in  three-dimensions,  the  cell  area)  and  its  location 
(centroid)  on  each  tilt.  The  data  obtained  for  the  lowest  tilt  is  assumed 
to  extend  to  the  surface.  The  height  of  the  centroid  is  determined  from 
the  range  to  the  centroid  and  the  elevation  angle  for  the  scan. 

The  three-dimensional  cell  or  volume  cell  is  defined  by  t lie  volume 
enclosed  within  the  peak  detected  cells  that  are  associated  from  scan-to- 
sean.  A criterion  similar  to  the  one  used  to  define  a peak  referenced 
cell  is  used  to  define  the  volume  cell,  from  the  sequence  of  attributes 
for  cells  detected  scan-to-scan , cells  1,  - and  3 on  this  figure,  a peak 
value  and  its  height  may  be  selected.  The  top  and  base  of  the  volume 
cell  are  defined  by  the  height  at  which  the  reported  average  value  drops 
the  required  number  of  quantization  units  below  the  peak  value.  Ibis 
height  is  determined  by  extrapolation  if  detected  cells  are  obtained  for 
each  tilt  or  is  taken  as  halfway  between  the  height  of  a value  above  the 
required  threshold  and  the  expected  height  of  the  intersection  of  the  cell 
with  the  next  tilt  plane  if  no  associated  cell  was  detected  on  the  next 
scan.  The  locations  of  the  peaks,  tops,  and  bases  are  depicted  on  the 
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figure.  The  volume  of  the  detected  volume  cell  is  assumed  to  be  the  pro- 
duct of  the  area  at  the  peak  value  (averaged  if  have  more  than  one  value 
equal  to  the  peak  value)  and  the  height  difference  between  top  and  base. 

The  volume  cells  may  have  a base  at  the  surface  as  illustrated  by 
cell  number  1 or  a base  above  the  surface.  In  the  latter  case,  the  cell 
may  be  detected  at  the  surface  as  illustrated  by  cell  7>  or  may  not  extend 
down  to  the  surface  as  illustrated  by  cell  2.  The  volume  cell  ts  des- 
cribed by  a set  of  attributes  which  includes  its  volume,  its  reflectivity 
and  area  at  the  surface  (zeroes  if  it  does  not  extend  to  the  surface) 
and  a number  of  additional  parameters  as  listed  in  Table  1.  In  some 
cases,  the  volume  cell  may  be  detected  on  some  scans  and  missed  on  inter- 
vening scans.  In  this  case,  the  volume  cell  will  be  "filled  in"  by 
interpolation  between  the  measured  values.  In  each  case,  the  centroid 
is  taken  to  be  the  surface  distance  at  the  lowest  elevation  angle  with 
a detected  cell  (single  scan). 

bach  detected  peak  is  tagged  with  the  identifier  of  the  enclosing 
lowest  threshold  fixed  contour.  The  fixed  contours  are  associated  from 
one  scan  to  the  next  using  the  identifiers  for  associated  peak  cells. 
Three-dimensional  fixed  contour  attributes  are  generated  using  the 
associated  data.  The  echo  tops  are  defined  as  lying  halfway  between  the 
height  of  the  highest  detected  peak  enclosed  by  the  contour  (which  may 
be  at  a threshold  level  many  dB  below  the  cell  top  reflectivity  level) 
and  the  height  of  the  cell  location  at  the  next  highest  elevation  angle. 
The  volume  within  the  fixed  contour  is  obtained  by  summing  the  areas  on 
each  scan  and  using  the  height  differences  between  the  centroid  locations. 
A complete  list  of  fixed  contour  volume  cell  attributes  is  listed  in 
Table  1. 

A. 2 Detection  and  Tracking  Algorithms 

The  volume  cell  attributes  are  obtained  from  the  peak  detected  cell 
attributes  obtained  on  each  azimuth  scan.  These  cells  must  be  associated 
from  scan-to-scan  before  the  data  from  different  scans  can  be  combined. 
This  association  process  lies  at  the  heart  of  the  volume  cell  detection 
and  tracking  process.  Association  is  affected  bv  pairing  cells  that  are 
sufficiently  close  together  to  come  from  the  same  three-dimensional  st rue 
ttire.  The  location  of  the  contours  surrounding  each  peak  is  not  available 


tor  uso  by  the  association  algorithm.  Association  is  accomplished  by 
first  selecting  cells  from  each  scan  whose  centroid  locations,  when  pro- 
jected on  the  surface,  are  separated  by  less  than  the  square  root  of 
the  combined  areas,  plus  the  distance  the  cell  might  move  during  the 
interval  between  scans,  plus  a fixed  distance  to  account  for  the  statis- 
tical centroid  location  uncertainty  introduced  by  the  expected  variability 
of  the  radar  data.  Comparison  is  always  made  between  the  cell  location 
for  a particular  scan  and  the  detected  volume  cell  location  from  the 
lowest  elevation  scan  on  which  it  was  observed. 

In  most  cases,  only  pairs  of  cells  will  be  evident  from  the  associ- 
ation algorithm  listed  above.  Occasionally,  in  a cluster  of  cells,  more 
than  one  cell  may  meet  the  association  criterion.  In  this  case,  the 
cells  will  be  associated  by  picking  those  cells  which  are  closest  after 
a uniform  offset  is  made  between  cell  locations  to  provide  the  largest 
number  of  associations.  This  is  illustrated  in  figure  3.  In  this  fig- 
ure, the  cell  to  be  tested,  cell  A,  is  closest  to  the  previously  detec- 
ted volume  cell  number  J,  but  after  translation  of  the  surface  locations 
relative  to  each  other  by  the  amount  shown  by  the  straight  line  on  the 
figure,  all  three  cells  (lettered)  may  be  associated  with  the  detected 
cells  (numbered)  with  smaller  physical  separations.  This  procedure  for 
conflict  resolution  (correlation)  will  be  used  only  when  the  resultant 
translation  is  physically  realisable  (i.e.  corresponds  to  a possible  cell 
tilt  or  translation  velocity). 

The  scan-to-scan  cell  association  logic  is  also  used  for  volume-to- 
volume  scan  sequence  association  for  cell  tracking.  Tor  the  tracking 
problem,  the  translation  analysis  for  conflict  resolution  becomes  more 
important  and  guidance  in  selecting  the  appropriate  translations  is 
taken  either  from  a prescribed  wind  vector  COO  mb  wind  for  instance)  or 
from  the  results  from  the  previous  scan.  1'he  conflict  resolution  or 
correlation  algorithm  is  performed  separately  for  large  echo  regions  in 
that  the  translation  values  are  allowed  to  vary  from  one  echo  region  to 
the  next. 

3.3  Volume  fell  \t tributes 

lhe  volume  cell  attributes  were  partially  defined  in  Section  3.1  and 
are  listed  in  fable  1.  \ different  set  of  attributes  is  generated  for 
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each  volume  cell  type.  Fixed  contour  attributes  include  measurements  of 
the  total  precipitation  produced  within  the  echo  region,  average  rain 
rates  on  the  surface,  profiles  of  area  and  of  average  reflectivity  for 
the  data  contained  within  the  fixed  contour  envelope,  and  summaries  of 
the  numbers  and  locations  of  peak  referenced  cells  detected  within  the 
echo  envelope.  The  peak  referenced  volume  cell  attributes  include  the 
size,  intensity,  and  location  parameters  listed  above  plus  information 
on  associated  volume  cells  of  different  types.  The  reflectivity  peak 
attributes  also  include  average  spread  and  shear  values  for  data  obtained 
within  the  cell  (between  base  and  top).  The  vertical  shear  values  are 
calculated  using  the  variation  of  the  average  radial  velocity  with  height. 

Fach  detected  volume  cell  is  characterized  by  its  attributes,  Each 
volume  cell  track  by  the  time  variation  of  its  attributes.  In  addition, 
tracks  may  be  characterized  by  a lifetime  and  average  velocity.  Obser- 
vations obtained  to  date  indicate  that  the  individual  peak  referenced 
reflectivity  cells  neither  merge  nor  split  but  only  develop  and  translate. 
Observations  of  a much  larger  data  set  are  needed  to  verify  this  model 
for  small  cell  behavior.  The  larger  fixed  contour  regions  may  merge  or 
split  depending  upon  its  stage  of  development.  By  tracking  the  enclosed 
cells  and  maintaining  the  cell  to  fixed  contour  identification  relation- 
ship, the  merger,  growth,  division  and  decay  of  the  larger  echo  regions 
may  be  automatically  taken  into  account.  The  morphology  of  the  echo 
development  process  may  be  important  to  the  understanding  of  precipita- 
tion dynamics,  however,  suitable  attributes  to  characterize  the  relevant 
processes  have  yet  to  be  defined.  Again,  experience  with  a significantly 
larger  data  set  is  required  before  further  progress  can  be  made. 

A. 4 Software  for  Volume  Cell  Detection  and  Tracking 

A computer  program  (TRACK)  was  generated  under  this  contract  to 
associate  cells  from  scan-to-scan,  resolve  conflicts  between  multiple 
associations  and  generate  the  volume  cell  attributes  listed  in  Table  I. 

The  listing  is  provided  in  Appendix  D and  operating  instructions  are 
given  in  Appendix  C.  The  program  processes  data  which  have  been  gener- 
ated by  the  cell  detection  program  (EXTRAD)  and  then  recorded  on  disk  or 
tape;  the  algorithms  described  in  Sections  3.1  to  3.3  are  used.  The  cur- 
rent output  is  a list  of  attributes  for  each  detected  volume  cell  (or 
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track).  These  outputs  are  again  stored  on  disk  or  tape  and  listed  on 
printout  for  subsequent  analysis.  Programs  have  not  been  generated  as 
vet  to  summarise  the  track  histories,  do  climatological  analysis,  or 
provide  samples  for  case  studv  analysis. 
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4.  CULL  UORUCASTI  Nil 


4.1  [Extrapolation  Along  Cell  Tracks 

Short  range,  0-20  minute,  forecasts  are  of  importance  for  severe 
jj  ' weather  warning  and  weather  hazard  avoidance.  The  objective  of  a fore- 

cast on  this  time  (and  by  implication  comparatively  small  distance)  scale 

I*  is  the  occurrence  of  a severe  event  at  a point  or  over  a small  area  such 

as  an  airport  or  the  approach  path  to  the  airport.  [Events  such  as  high 
winds  produced  by  downbursts  or  by  gust  fronts,  or  regions  of  intense 
hail  fall  are  of  most  interest.  Larger  time  and  area  forecasts  may  be 
made  of  the  probable  occurrence  of  a severe  weather  event  but  the  pro- 
blem considered  here  is  the  remote  detection,  tracking,  and  prediction 
of  the  development  and  decay  of  a severe  weather  event.  The  first  parts 
of  the  problem,  the  detection  and  tracking  of  small  scale  features  in 
the  radar  data,  were  considered  in  previous  sections.  The  association 
of  these  features  with  severe  weather  events  is  a second  problem  that  is 
not  being  considered  under  this  contract.  The  third  part,  the  forecast 
of  cell  development  and  motion  are  given  initial  consideration  in  this 
sect  ion . 

The  short  range  forecast  problem  considered  here  is  different  from 
most  meteorological  forecast  problems  in  that  the  occurrence  of  at  least 
one  cell  is  assumed.  The  problems  are  where  will  it  move,  how  will  it 
develop  in  time,  where  will  new  ones  develop,  and  will  it  produce  a severe 
weather  event.  The  measure  of  success  of  the  forecast  procedure  must 

I depend  upon  the  particular  problem  addressed.  The  success  of  a cell  posi- 

tion forecast  should  he  measured  by  the  distance  between  its  actual  posi- 
tion and  forecast  position  at  forecast  time.  The  success  of  the  cell 
development  forecast  should  he  measured  by  the  difference  between  the 
actual  intensity  (reflectivity,  say)  and  forecast  intensity  at  forecast 
time.  Classical  evaluation  procedures  that  test  the  occurrence  or  non- 
• occurrence  of  the  event  at  a number  of  geographical  locations  are  not 

recommended  in  that  the  reason  for  failure  of  the  forecast  may  not  be 
readily  apparent. 

Crane  (197b)  used  data  obtained  by  AIM.  at  Wallops  Island  to  perform 
some  initial  tests  of  cell  tracking/track  position  forecast  algorithms, 
tie  assumed  that  a single  cell  motion  vector  would  adequately  describe  the 
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propagation  of  a field  of  cells  and  tested  that  assumption  hy  measuring 
the  along  track  and  cross  track  position  errors  as  a function  of  forecast 
time  (0-20  minutes).  His  results  showed  that  the  forecast  error  (posi- 
tion) was  the  order  of  the  cell  diameter  (median  value)  at  the  half  life 
(median  lifetime)  of  the  cell.  He  also  found  that  the  half  life  was  of  the 
order  of  10  to  15  minutes  and  the  median  cell  diameter  was  of  the  order  of 
5 km.  This  result  was,  however,  based  on  a limited  amount  of  data. 

Crane  found  that  the  operation  of  the  track  extrapolation  forecast 
procedure  could  be  improved  by  using  different  motion  vectors  for  differ- 
ent regions  of  the  display  area  (see  figure  4 for  an  example  showing  the 
different  directions  possible  for  cell  tracks  observed  during  the  same 
time  interval).  This  idea  is  incorporated  in  the  volume  cell  detection 
and  tracking  program.  The  major  problem  of  using  extrapolation  along  a 
track  for  forecasting  is  the  relatively  short  lifetime  of  most  cells. 

The  largest  and  most  intense  cells  persisted  for  a relatively  long  time, 

40  to  50  minutes.  Except  for  these  cells,  the  cell  could  disappear  by 
forecast  time  and  a significant  number  of  new  cells  could  appear,  fore- 
casts based  upon  cell  time  histories  could  be  used  to  estimate  cell  life- 
time but  the  forecast  of  new  cell  site  development  is  a major  problem  for 
any  radar  data  based  short  range  forecast  procedure. 

4.2  New  Cell  Site  Development 

The  requirement  for  a procedure  to  forecast  the  most  probable  sites 
for  new  cell  development  and  the  observation  of  a seemingly  persistant 
organization  or  structure  for  the  location  of  active  cells  lead  to  an 
investigation  of  the  most  probable  distance  or  spacing  between  the  cells. 
The  Wallops  Island  data  previously  processed  by  Crane  (1‘C'o)  were  used  for 
this  analysis.  The  data  consisted  of  computer  prepared  maps  of  detected 
cell  locations  for  each  of  the  azimuth  scans  processed  from  the  ID~5 
summer  data  set  provided  by  API..  Nearest  neighbor  distances  were  obtained 
from  the  data  to  investigate  the  existence  of  persistant  dominant  scales 
for  cell  organization.  The  smallest,  second  and  third  shortest  distances 
between  each  cell  and  its  neighbors  were  determined  and  tabulated  for 
each  cell  for  each  available  volume  scan.  It  is  noted  that  with  the  pro- 
cessing technique  used,  the  distance  between  two  closely  spaced  cells  is 
counted  twice,  each  cell  being  the  nearest  neighbor  of  the  other.  Nearest 
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neighbor  distributions  were  generated  for  each  scan.  The  mode  for  each 
distribution  was  used  to  estimate  the  most  probable  cell  spacing  for  that 
scan.  Time  histories  of  the  modal  values  for  the  smallest,  second  and 
third  shortest  distances  are  given  in  Figure  5 for  three  days  with 
relatively  large  numbers  of  cells  (more  than  20  per  scan).  These  data 
show  nearest  neighbor  distances  to  be  on  the  order  of  8 km.  Some  perio- 
dicity is  evident  in  the  nearest  neighbor  distance  values  but  this  is 
most  likely  due  to  the  small  number  of  samples  used  to  generate  a dis- 
tance estimate.  The  data  show  a marked  day  to  day  consistency. 

The  most  likely  sites  for  the  development  of  new  cells  in  a cluster 
of  cells  should  be  at  the  nearest  neighbor  spacing  from  the  extrapolated 
position  of  observed  cells.  A preliminary  test  of  this  hypothesis  is 
illustrated  in  Figures  6 and  7.  Figure  6 displays  the  locations  of  new, 
short-lived  or  decaying,  and  trackable  cells  observed  at  2229  GMT  on  18 
June  1973  at  Wallops  Island.  The  trackable  cells  are  depicted  at  mid- 
track. Since  the  observations  are  separated  by  more  than  20  minutes, 
only  long-lived  cells  are  trackable  and  the  majority  of  cells  are  depic- 
ted as  new  (first  detection  at  2229  GMT)  or  short-lived  (detected  on  the 
previous  scan).  Figure  7 displays  the  situation  more  than  an  hour  later, 
at  2337  GMT.  The  new  cells  depicted  at  this  time  are  designated  by  N. 

The  older  cells  shown  on  Figure  6 have  been  extrapolated  to  the  time  for 
Figure  7 using  the  mean  cell  velocity  obtained  from  the  trackable  cells. 

It  is  noted  that  most  of  the  cells  have  decayed  by  this  time  and  the 
locations  represent  the  extrapolated  locations  of  the  cell  sites.  The 
circles  drawn  about  each  of  the  extrapolated  cell  sites  have  radii  equal 
to  the  most  likely  nearest  neighbor  distance.  The  results  show  that  23 
(66°s)  of  the  new  cells  appear  at  or  within  2 km  of  the  nearest  neighbor 
distance  while  12  lie  outside  this  range.  These  results  show  consider- 
able promise  given  the  relatively  long  time  intervals  between  observations. 

The  new  cell  sites  depicted  on  Figure  7 also  show  a number  of  cells 
approximately  spaced  from  each  other  by  the  nearest  neighbor  distance  but 
not  connected  to  an  older  cell  site.  A forecast  based  upon  a more  com- 
plex pattern  vising  the  nearest  neighbor  distances  in  a more  regular  ex- 
tended structure  may  be  capable  of  predicting  the  locations  of  a larger 
number  of  cells.  The  structure  may  also  reduce  the  size  of  the  region 
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forecast  for  probable  new  cell  development  by  specifying  a segment  or 
segments  of  the  nearest  neighbor  arc.  A considerably  larger  s;uiiple  of 
data  with  smaller  time  steps  between  observations  is  required  before  a 
procedure  for  new  cell  site  forecasting  can  be  developed  and  tested. 

The  preliminary  results  presented  here  are  encouraging  and  indicate  that 
a new  cell  site  forecasting  procedure  can  be  developed. 


5.  ANALYSIS  OF  SELECTED  DATA  SETS 


5.1  Use  of  the  Computer  Programs 

The  computer  programs  described  above  were  developed  as  a prototype 
for  a real-time  prc  sing  system  to  be  employed  in  reducing  the  number 
of  data  bits  required  to  describe  a set  of  radar  observations.  The  ob- 
jective of  the  program  was  to  replace  the  large  volumes  of  data  obtained 
by  a radar  system  to  a set  of  fixed  contour  and  peak  referenced  cell 
attributes  capable  of  representing  the  same  essential  information.  A 
prototype  program  is  now  available  with  optional  choices  for  many  of  the 
processing  parameters  and  for  representing  the  cell  and  contour  attributes. 
The  major  thrust  of  this  contract  has  been  program  development  and  ini- 
tial testing  of  critical  hypotheses  such  as  the  possibility  of  forecast- 
ing the  locations  of  new  cell  sites.  A considerable  amount  of  work  still 
remains  using  the  programs  to  optimize  the  processing  parameters. 

An  example  of  the  input  radar  data  for  a pair  of  fixed  contours  is 
given  in  Figure  1.  The  raw  data  are  the  three  input  parameters,  reflec- 
tivity, radial  velocity,  and  Doppler  (or  velocity)  spread.  An  internally 
generated  data  field,  tangential  shear,  is  also  depicted.  Contoured  data 
for  the  three  fields  used  in  further  processing  for  one  of  the  echo  regions 
are  depicted  in  Figure  8.  In  this  figure,  the  quantization  steps  for 
each  data  field  are  1 dB,  0.5  m/s,  and  0.1  m/s/km  for  the  reflectivity, 
spread,  and  tangential  shear  fields  respectively.  The  reflectivity  data 
show  relatively  long,  thin  contours.  The  basic  resolution  element  for 
processing,  1°  by  0.5  km,  is  a square  or.  this  figure.  For  more  efficient 
processing,  the  resolution  elements  should  be  adjusted  so  the  contours 
to  be  detected  are  nearly  circular.  The  radar  data  are  read  into  the  com- 
puter in  1°  by  0.15  km  resolution. elements.  In  this  example,  the  data 
along  a radial  should  be  averaged  by  at  least  a factor  of  6 (5  times  the 
factor  used  for  this  figure)  to  produce  1°  x 0.9  km  resolution  elements. 
This  adjustment  may  be  made  automatically  in  the  cell  detection  program. 

The  spread  and  tangential  shear  data  show  large  changes  from  one  resolu- 
tion element  to  the  next  although,  with  the  quantization  steps  used,  the 
data  reveal  the  same  tendency  toward  elongation  as  the  reflectivity  data. 

Contours  enclosing  peak  detected  cells  are  indicated  by  the  cross 
hatching  on  Figure  8 and  on  Figure  9.  The  contours  are  generated  at 


The  out  or  oontouv  for 


throo  quantization  units  below  t ho  peak  values, 
reflectivity  on  figure  S arul  the  outer  contour  on  figure  9 correspond  to 
the  JO  dBZ  threshold  used  for  processing.  Two  reflectivity  peaks  were 
detected  although  a shoulder  on  the  first  peak  (in  range,  the  shoulder 
is  at  shorter  ranges)  would  have  been  detected  as  an  individual  cell  if 
two  quantisation  steps  were  used  to  define  the  cell.  1'he  shoulder  would 
also  he  included  with  the  two  detected  cells  in  a single  cell  if  the  1" 

In  0,9  km  resolution  area  were  used. 

The  tangential  shear  and  Poppler  spread  peaks  are  shown  superimposed 
on  the  reflectivity  contours  in  figure  9,  The  two  tangential  shear  peaks 
occur  for  each  reflectivity  peak  if  the  shoulder  were  included.  These 
peaks  together  with  the  reflectivity  data  appear  to  be  delineating  the 
updraft,  down  draft  structure  of  the  cell.  Kith  a change  to  a 1''  hv  0,9 
km  resolution  area,  these  sub-cell  structures  would  be  averaged  to  a con- 
sistent, one  cell  picture  for  this  echo  region.  The  Poppler  spread  data 
reveal  two  types  of  peaks,  four  of  the  seven  detected  peaks  coincide  with 
the  tangential  shear  peaks  outlining  regions  of  high  shear  and  high  spread. 
These  peals  are  also  intimately  related  to  the  reflectivity  structure. 

Three  other  Poppler  spread  peaks  also  are  evident  at  the  von  edges  of 
the  JO  dBZ  contour.  These  peaks  correspond  to  regions  of  relatively  low 
reflect  ivit\  (..'a  dBZ  and  lower)  and  austonsiblx  represent  regions  where 
low  s ignal - to-no i se  values  give  rise  to  deceptivel\  high  Poppler  spreads. 

The  number  and  shapes  of  the  peak  referenced  cells  will  var\  as  the 
quantization  step  size,  the  number  of  quantization  steps  for  peak  detec 
tion,  anil  the  si.-e  of  the  resolution  area  are  changed.  More  experience 
is  required  using  these  computer  programs  to  select  the  optimum  contbina 
tion  of  these  parameters  for  the  detection  of  phvsic.illx  meaningful  cells. 
As  used  in  generating  figures  S and  9,  it  appears  that  too  much  detail 
is  present  and  some  of  the  detected  regions  especially  for  spread  and 
tangential  shear  data  represent  structure  on  too  small  a scale  (for 
example,  the  individual  updraft  and  downdraft  regions  within  a cell). 

5.J  Volume  fell  Observations 

The  volume  cell  detection  program  is  used  to  combine  the  data  from 
a number  of  individual  azimuth  scans.  Pata  obtained  from  the  Sudbur> 
radar,  subsequent  to  19JS  OMT  on  August  lo,  19**S  using  1''  In  ('.9  km 


resolution  areas  are  depicted  in  Figures  10-14.  These  data  show  the  cen- 
troid location  of  each  reflectivity  peak  referenced  cell  coded  both  by 
average  reflectivity  and  by  area.  The  locations  of  the  tangential  shear 
peaks  are  denoted  by  X and  the  locations  of  Doppler  spread  peaks  are 
denoted  by  +. 

The  data  from  each  of  the  scans  were  combined  to  form  volume  cells 
shown  in  Figure  15.  In  this  figure,  the  solid  lines  connect  the  cell 
locations  as  detected  on  each  scan  using  the  algorithms  described  in 
SEction  3 (Track  Program).  Cells  within  50  km  of  the  radar  were  detected 
at  elevation  angles  up  to  5°;  cells  at  further  ranges  were  only  detected 
at  elevation  angles  below  4°  (heights  less  than  5 km).  A limited  set  of 
volume  cell  attributes  for  this  data  set  is  listed  in  Tahle  2. 


I inure  11  Reflectivity,  Tangential  Shear  (X),  and  Doppler  Spread 
Peak  Locations  for  a 2°  ! location  Angle 
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Reflectivity,  'tangential  Shear  fXj,  and  Lioppler  Spread 
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6.  Rl-COMMliNUAT  IONS 

6.1  Parameter  Optimization 

A set  of  computer  programs  has  been  generated  to  provide  an  auto- 
matic means  for  the  extraction  of  information  from  large  volumes  of  radar 
data.  The  programs  are  written  to  be  as  general  as  possible  to  enable 
rapid  changes  in  processing  parameters.  As  indicated  in  Section  5,  the 
optimum  values  for  these  parameters  are  not  known  and  must  be  determined. 

It  is  expected  that  the  parameters  should  change  from  one  radar  system 
to  the  next  depending  principally  on  the  resolution  volume  and  number  of 
independent  samples  per  resolution  element.  The  first  problem  to  be  con- 
sidered in  the  use  of  this  set  of  programs  is  the  optimization  of  para- 
meters. This  can  only  be  done  by  processing  a relatively  large  number  of 
radar  scans  for  different  rain  conditions. 

Ideally,  auxiliary  data  should  be  available  to  provide  a standard 
of  comparison  for  the  output  from  the  program.  Many  case  studies  such 
as  the  one  reported  in  Section  5.1  should  be  performed  to  obtain  the 
raw  contour  data  to  provide  a comparison  standard.  For  use  in  detect- 
ing severe  weather  events,  auxiliary  data  on  the  severe  events  are  also 
requ i red . 

6.2  Real  lime  Processing 

The  programs,  although  general  in  nature,  were  written  with  the 
ultimate  goal  of  use  in  a real  time  processor.  After  the  analysis  for 
parameter  optimization  has  been  completed,  specialization  to  a real  time 
processor  may  be  accomplished.  Real-time  processing  requires  the  minimi- 
zation of  computer  storage  and  operating  time.  Major  steps  can  be  made 
in  this  direction  by  increasing  the  resolution  element  area  for  process- 
ing (averaging  over  M km  in  range  as  recommended  in  Section  5.1)  and  bv 
reducing  the  number  of  peak  detection  operations.  The  latter  can  be 
accomplished  by  processing  either  tangential  shear  or  Popp l or  spread  data 
but  not  both  as  is  currently  done.  Another  time  saving  step  is  to  reduce 
the  volume  of  output  by  not  preparing  the  fixed  level  contour  plotting 
displays  in  the  computer  but  doing  the  fixed  contouring  in  the  color  dis- 
play processor  as  is  currently  done  for  the  output  of  digital  integrators. 
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Preprocessing  of  the  data  for  calibration,  velocity  ambiguity  resolution 
and  conversion  from  variance  to  velocity  spread  will  also  save  some  time 

Current  running  time  for  reflectivity  processing  only  but  including 
fixed  contour  generation  for  display  is  2 minutes  per  azimuth  scan.  By 
just  doing  the  preprocessing,  the  running  time  could  be  reduced  by  more 
than  a factor  of  two.  The  other  reductions  recommended  above  plus  inter 
nal  programming  changes  to  reduce  the  use  of  indirect  array  addressing 
should  result  in  a program  that  will  do  at  least  two  data  fields,  reflee 
tivity  and  tangential  shear  or  Doppler  spread,  in  real  time  (30  seconds 
per  scan) . 

(i.3  Spatial  Analysis  of  fell  Development 

The  initial  considerations  of  the  forecast  of  new  cell  site  loca- 
tions indicates  that  improvements  should  be  possible  if  attention  is 
focused  on  the  structure  or  organization  of  the  cell  location  patterns. 
Processing  to  date  has  used  the  determination  of  nearest  neighbor  dis- 
tances to  obtain  information  about  cell  spacings.  The  nearest  neighbor 
distances  provide  estimates  of  the  locations  of  secondary  maxima  in  the 
spatial  correlation  function  for  cell  locations.  Information  on  struc- 
ture can  better  be  obtained  from  more  refined  correlation  function  (or 
spatial  power  spectral  of  cell  location.  These  analyses  should  be  con- 
ducted using  a much  larger  data  sample  preferably  for  a number  of  differ 
ent  locations  and  storm  types. 

6.4  Morphological  and  Climatological  Analysis 

The  programs  described  above  are  used  to  extract  the  significant 
information  from  a large  volume  of  radar  data.  The  result  is  still  a 
formidable  data  set  comprised  of  a number  of  lists  of  volume  cell  attri- 
butes for  each  scan,  storm,  and  day.  These  data  must  in  turn  be  reduced 
to  a manageable  set  to  describe  the  morphology  of  cell  development.  Thi 
analysis  must  be  performed  before  any  meaningful  cell  forecast  procedure 
can  be  developed  and  tested.  The  analysis  entails  both  the  construction 
of  new  programs  and  the  processing  of  large  volumes  of  data. 
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APPENDIX  A 


CELL  DETECTION  AND  TRACKING  PROGRAM 
INSTRUCTIONS  FOR  OPERATION 

A.l  Description  of  Input  and  Output 

Program  input  and  output  are  depicted  in  Figure  AI . The  tape  input 
format  is  given  in  Table  Al.  The  control  cards  are  discussed  in  section 
A2.  The  program  produces  (a)  tapes  of  computed  attributes  for  input  to  a 
second  program  for  computing  volume  scans;  (b)  a plot  tape  is  generated 
that  can  be  stored  for  input  to  another  program  "EXPAND"  which  is  a gen- 
eral purpose  plotting  package  for  plotting  the  fixed  contours,  centroids, 
cell  identification  and  peak  locations  expanded  over  selected  areas;  (c) 
B-scan  maps  are  also  produced  as  an  option;  and  (d)  at  the  completion  of 
a scan  the  program  will  print  out  fixed  contour  attributes,  peak  detected 
cell  attributes  and  tangential  shear  maxima  attributes.  All  of  the  attri- 
butes printed  have  identifiers  which  can  be  associated  with  the  identi- 
fiers displayed  on  the  expanded  plots. 

A. 2 Control  Card  Format 

Control  card  input  to  the  program  is  NAMELIST  input  which  allows 
certain  parameters  in  the  program  to  default  or  to  be  set  to  different 
values.  The  variable  names,  type  (LOGICAL  L,  INTEGER  I,  and  REAL  R) , 
dimension,  default  value  and  their  meanings  are  listed  in  Table  A2. 


Dopple 


12  Bit  Wo 
Position 


Hf°°- 


Least  Significant  Bit 


tabu:  a: 


CARO  FORMAT  FOR  PROGRAM  EXTRAD 
Reads  in  program  parameters  via  NAMELIST  format. 


NAMELIST  VARIABLES:  (.Level  7804151 


NAME 

TYPE. 

01  MENS  I ON 

DEFAULT 

MEANING 

PR INTI 

L 

1 

FALSE 

Unused . 

print: 

L 

1 

FA  LSI- 

If  .TRUE,  print  B-scan  maps  of 
mean  and  variance. 

PR  I NTS 

L 

1 

FALSE 

If  .TRUE,  print  B-scan  maps  of 
dB;. 

PRINT4 

L 

1 

FALSE 

Unused . 

I COPES 

I 

3 b 

Blank  thru 

Z then  1 
thru  9. 

Codes  for  representing  dB;  cate- 
gories for  B-scan  map  output. 

Al 

R 

1 

. 13779 

In  the  linear  equation  y = mx+b 
for  computing  coded  dBz  for  B- 

R1 

R 

1 

1.5 

scans,  Al  = M and  B1  = b. 

A2 

S 

V 

1 

.017 

Unused . 

B2 

R 

1 

1 8 . () 

Unused . 

CONTRZ 

L 

1 

TRUE 

If  FALSE,  do  not  process  peak 
cell  attributes. 

CONTRV 

L 

1 

TRUE 

If  FALSE,  do  not  process  reflec- 
tivity, shear,  ti  spread. 

CONTRS 

L 

l 

TRUE 

If  FALSE,  do  not  process  shear 
and  spread. 

NFII.E 

I 

1 

0 

Number  of  files  on  tape  to  skip 
before  processing. 

NOME 

I 

1 

1 

Number  of  files  on  tape  to  process 

AC 

R 

4 

-107.7, 

♦1.97, 

-0.94, 

+.0018 

Calibration  coefficients  for 
computing  DBM  below  a threshold 
XCIIT.  (see  XCUT) 

NAME 

COPLOT 

VOLTAP 

CALM 

CALB 

XCUT 

CK 

ZMAX 

VMAX 

NREC 

NIIMR 

IRON 

INC 

STARTR 

STOPR 

INPRF 

SCALE 

AH 

AA 

BB 


TYPE 

DIMENSION 

DEFAULT 

MEANING 

L 

1 

FALSE 

If  TRUE,  output  a tape  for 
plotting. 

L 

1 

TRUE 

If  FALSE,  do  not  output  an 
attribute  tape. 

R 

1 

. 332 

In  the  calibration  equation  for 

R 

1 

-98.3 

DBM  v = mx+b,  CALM  = M t)  CALB  = b. 

R 

1 

10.0 

Threshold  value  that  determines 
which  equation  to  use  for  cali- 
bration. (linear  or  non-linear) 

R 

1 

10.0 

In  the  equation  for  computing  dB; , 
K+P+20ALOG 10(S(I-.S))  K = CK. 

R 

1 

0.0 

Not  currently  used. 

R 

1 

0.0 

Not  eu r ren 1 1 y us ed . 

I 

1 

1 

Not  currently  used. 

I 

1 

999 

Number  of  radials  to  be  processed. 
Use,  default  value  when  doing  full 
scan. 

I 

1 

0 

Run  number  chosen  by  user . 

I 

1 

0 

Not  currently  used. 

R 

1 

0.0 

Where  along  a radial  in  kilometers 
processing  is  to  start. 

R 

1 

300.0 

Where  along  a radial  in  kilometers 
processing  is  to  stop. 

I 

1 

3333 

Value  of  PRF  (Pulse  Repetition 
Frequency)  to  be  used  when  PRF 
cannot  be  obtained  from  the  data 
tape. 

R 

I 

1.0 

Scale  factor  for  drawing  fixed 
contours . 

R 

1 

1.21 

Constant  for  computing  heights 
of  cells. 

R 

1 

300 

Constant  for  computing  rain  rate. 

R 

1 

1.5 

Exponent  for  computing  rain  rate. 
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NAME 

'IT  PL 

DIMENSION 

DEFAULT 

MEANING 

XI 

R 

1 

0.0 

Frame  size  coordinates  for  fixed 
contour  plotting.  Less  than  or 
equal  to  8 inches. 

X2 

R 

1 

8.0 

Same  as  above. 

Y1 

R 

1 

0.0 

Same  as  above. 

Y2 

R 

1 

8.0 

Same  as  above. 

TV 

T 

1 

3S 

Mean  wind  velocity  in  a fixed 
echo  contour  is  not  computed  for 
dBz  greater  than  this  value. 

TSV 

R 

1 

iob 

Not  currently  used. 

L.PV 

I 

1 

3 

Cell  detection  threshold  for 
reflectance  peaks. 

LTV' 

I 

1 

3 

Cell  detection  threshold  for 
velocity  peaks. 

LSV 

I 

1 

3 

Cell  detection  threshold  for 
shear  peaks. 

ICOMP 

I 

1 

6 

Data  compression  factor,  (range 
integration! 

VMIN 

R 

1 

0.0 

Unused. 

SVMIN 

R 

1 

0.0 

Unused . 

ARF.AMN 

R 

1 

1.0 

Any  completed  contour  having 
an  area  less  than  AREAMN  will 
be  ignored  but  is  plotted  if  a 
fixed  contour. 

WAV  EL 

R 

1 

0.0542 

Radar  wavelength  in  meters. 

VQUANT 

R 

l 

10.0 

Tangential  shear  quantization 
step  ( m/s/km !'  . 

SQIIANT 

R 

1 

2.0 

Doppler  spread  quantization 
step  (m/s)  . 

RQIIANT 

R 

1 

1.0 

Reflectivity  quantization  steps 
(indB-1). 

DAZM 

R 

1 

1.0 

Beam  width  (degrees). 

f 
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p*  - ' 

1 

... 



1 

NAME 

TYPE 

DIMENSION 

DEFAULT 

MEAN  1 NG 

ESTART 

R 

1 

(1.5°) 

Elevat ion 

start  angle 

(degrees) . 

DELT 

R 

1 

(.5°) 

Delta  elevation  angle 
defining  next  scan. 

(degrees) 

L 
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APPENDIX  B 

PLOTTING  PROGRAM  'EXPAND'  (VERSION  1.0) 

B.l  Description  of  Input  and  Output  /' 

/ 

Program  EXPAND  utilizes  the  plot  tape  generated  by  program  'EXTRAD' 
as  input  to  generate  not  only  full  scan  plots  of  fixed  contours  and 
their  centroids  but  on  option  will  expand  and  plot  certain  areas  of 
interest.  Also  on  option,  it  will  plot  out  locations  of  centroids  of 
fixed  contours,  peak  detected  cells,  Doppler  spread  and  tangential  shear 
locals.  These  plotting  options  apply  to  full  scan  plots  as  well  as 
expanded  plots. 

The  plots  are  generated  on  the  CALCOMP  ink  pen  drum  plotter.  The 
X axis  are  labeled  negative  kilometers  west  of  the  radar  and  positive 
east.  Y axis  are  labeled  negative  kilometers  south  of  the  radar  and 
positive  north.  The  date  and  elevation  angle  are  also  annotated. 

B.2  Control  Card  Format 

Control  card  input  to  the  program  is  NAMELIST  input  which  allows 
certain  parameters  to  default  or  to  be  set  to  different  values.  The 
variable  names,  type,  dimension,  default  and  meaning  are  listed  in  Table 
Bl. 

It  is  suggested  that  program  EXPAND  generate  full  scan  plots  of 
fixed  contours  from  the  entire  EXTRAD  tape  first  before  generating 
expanded  areas  of  view.  This  allows  one  to  examine  exactly  what  each 
scan  contains  and  where  expansion  would  be  of  interest.  One  set  of 
NAMELIST  input  cards  are  needed  for  each  scan.  If,  for  example,  it  is 
desired  to  go  into  the  third  scan  on  tape,  three  NAMELIST  set  ups  must 
oc  cur . 
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CARP  FORMAT  FOR  PROORAM  EXPAND 
READS  IN  PROORAM  PARAMETERS  VIA  NAMEl.IST  INPUT 


N AMP. 1. 1ST  VAR l AHI, I S: 


NAMP 

PYI’P 

DIMENSION 

mi  Aiu/r 

MEAN  INC 

[I'l.T 

4 

m.T(i) 

. PAl.SP. . 

If  .TRIIP.  plot  a Jot  to  locate 
peak  detec t ed  cel  Is. 

IPLT(2) 

.FA  l.S  I:. 

If  .TIUIP..  plot  a C to  locate 
fix  contour  centroids. 

1 IMjT(.V) 

. PALSP . 

If  .TRUli.  plot  a ♦ to  locate 
Doppler  spread. 

I PI.  1(1 ) 

. FALSE. 

If  .TIUIP.  plot  an  X to  locate 
tangential  shear. 

z\ 

1 

.TRUE . 

Plot  fix  contours  for  first 
level  contour. 

•y  -y 

1 

. TRIIP . 

Plot  fix  contours  for  second 
level  contour. 

23 

1 

.TRUE. 

Always  set  to  .PAl.SP.. 

24 

l 

.TRIIP. 

Mwavs  set  to  .PAl.SP,. 

l.S 

'■ 

I 

. P Al.SP . 

When  .TIUIP!.  expanded  area  plot 
are  requested.  When  .PAl.SP. 
draw  full  scan  only. 

XKI 

R 

1 

-.A  20.0 

Western  plot  area  limit  (km). 

xk: 

R 

1 

320.0 

Pastern  plot  area  limit  (km). 

YKl 

R 

1 

-320.0 

Southern  plot  area  limit  (km). 

YK2 

R 

1 

320.0 

Northern  plot  area  limit  (km). 

I.K 

•• 

1 

.PAl.SP. 

If  .TIUIP!.  plot  centroids  of 
given  areas. 

*Note : 

Expanded 

plots  will  always 

square  off 

any  rectangular  area  request 

to  the  largest  ;ixis  requested. 
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TRACK  I N('.  PROCRAM  IASOCCI.) 

C.l  I'escr  ipt  ion  of  Input  ami  Output 

I'he  input  tape  is  the  tape  produced  by  program  "IXTRAO".  This  pro 
gram  takes  the  roll  attributes  ami  associate  cells  for  tracking  volume 
cells.  The  output  listing  includes  the  locution  of  volume  cells  and  th 
updated  cell  II'. 

C.J  Control  Card  format 

Control  card  input  to  the  program  is  NAM1T.1ST  input  which  allows 
certain  parameters  in  the  program  to  default  or  to  be  set  to  different 
values.  The  variable  names,  type  (.I.OCICAI.  I.,  INTLC.LR  I,  and  RIAL  Rt, 
dimensions,  default  values  and  their  meanings  are  listed  in  fable  Cl. 


TABl.r  Cl 


CARP  IORMA  1'  I'OR  PROGRAM  ASOCCI 


Reads  in  Program  Parameters  via  N AMI  11  SI"  format 
NAM1  1 1ST  VAR  1 AIM. IS: 


N AMI 

TV  PI 

m Ml  NS  ION 

niiAiii.r 

Ml  \NINC. 

PR  l 

i 

1 

1 AI.SI 

When  .TRIM',  program  prints 
out  the  input  attributes  data. 

PR  J 

i 

1 

TAI  SI 

When  .'TRIM  . program  prints  out 
the  initial  volume  cell  attri 
luites  and  the  total  number  of 
cells. 

PR  A 

l 

l 

TAI.S1 

When  .TRIM'.  program  prints  out 
the  associated  update  cell 
attributes  and  the  total  nnm 
her  of  associated  cells. 

PR-1 

i 

l 

i'ALSI 

Current 1\  unused. 

im  vi 

i 

1 

1 

1 nput  t ape  or  f i 1 e . 

ini  v: 

i 

1 

1 

Current l \ unused . 

1 STOP 

i 

l 

1 

Stop  progi am,  when  "0"  continue 
process i ng . 
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